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Abstract 

In the mode coupling theory of the liquid to glass transition the long time 
structural relaxation follows from equations solely determined by equilibrium 
structural parameters. The present extension of these structural relaxation 
equations to arbitrarily short times on the one hand allows calculations un- 
affected by model assumptions about the microscopic dynamics and on the 
other hand supplies new starting points for analytical studies. As a first appli- 
cation, power-law like structural relaxation at a glass-transition singularity 
is explicitly proven for a special schematic MCT model. 
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§1. Introduction 



Dense liquids exhibit anomalously slow, temperature sensitive and non-exponentional 
relaxation processes. They are conventionally termed 'structural relaxation' because cooper- 
ative rearrangements of particles are presumably involved. The mode coupling theory (MCT) 
suggests that the evolution of the anomalous dynamics and its splitting off from the normal 
liquid dynamics is caused by bifurcation singularities in non-linear retarded equations of 
motion for the (normalized) intermediate scattering functions ^q{t) = Sq{t)/Sg (Leutheusser 
1984, Bengtzelius et al. 1984). These functions are some of the most simple ones specifying 
structural dynamics. The bifurcations of MCT, and the long-time dynamics in their vicinity, 
are solely determined by the equilibrium (static) structure factor Sq, see (Gotze and Sjogren 
1989, Pranosch et al. 1998), and thus the MCT provides a frame to single out and to define 
structural relaxation. In this contribution the MCT equations for the structural relaxation 
are extended to arbitrarily short times thus rendering them well defined even in the limit 
that the microscopic transient dynamics can be neglected. Whereas in (Pranosch et al. 1998) 
this was achieved by a discrete-dynamics model, here integral equations with explicit ini- 
tial conditions are formulated. Thus the dependence on the initial conditions is displayed 
clearly. Additionally, monotone solutions are guaranteed (Fuchs et al. 1991). Furthermore, 
for the obtained structural dynamics equations non-exponential relaxation can be proven 
for a special MCT model. 



The MCT equations of motion for the autocorrelation functions of density fluctuations in 
simple liquids, ^q{t), are derived from Newton's equations with the Zwanzig-Mori formalism 
leading to: 



§2. The equations of motion of MCT 



di%{t) + nl%{t) + /o* dt' Mq{t - 1') dt,%{t') = , 



(1) 



$,(i^0) = l-i(Q,i)2, 



(2) 
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and subsequent mode coupling approximations for the generalized longitudinal viscosity 
finally resulting in: 

M,(^) = ^]2^,(t) + M-g•(^), (3) 

nig{t) = Ek+p=q ^(q; k, p) Mt)Mt) ■ (4) 

In eqn. (|I]), the frequency Qq sets the microscopic time scale for the bare oscillators with 
wave vector q. In eqn. viscous damping arises from mode-mode coupling, mq(t), eqn. 
(I), and from microscopic short-time dynamics denoted by ikf'^S'; see (Gotze 1991) for 
explicit expressions and discussions. A Laplace transformation of eqs. (|l|) and (|) using 
<|)(s) = dte-'' ^t) leads to 

l>(s) = 1 / [ s + 1 / [ m{s) + is + W'^- is))/nl ]]. (5) 

The existence of a unique solution of the MCT equations which has all properties of an 
autocorrelation function has been shown by Haussmann (1990). Further and more detailed 
properties could be proven by Gotze and Sjogren (1995) upon the assumption of overdamped 
short time motion. Let us mention two results relevant for this work, namely (i) the existence 
of bifurcations in the equations derived from eqs. (^ and (^) for the long-time limits, fq = 
limt_,oo '^'i3(^); here the fg are the glass form factors and the bifurcations happen for critical 
values of the mode coupling vertices, V = V^, where fq = f^. {ii) for regular vertices, the fq 
depend smoothly on the V , and there exists a longest relaxation time, l/F, cutting off the 
final relaxation exponentially, $g(t) — fq = 0{e~^^) for V ^ V^. 

One of the most simple models exhibiting the generic features of the MCT bifurcation 
dynamics in eqs. (|l]) to (^) is obtained by the simplification to study only one correlator, $(t), 
which experiences non-linear feedback via one mode-coupling functional m(t) = fi$(t) + 
V2$^(t) with fi 2 > (Gotze 1984). We will start by discussing this so-called Fi2-model 
before generalizing our results to the equations relevant for simple liquids. 



3. Results for the schematic F12— model 
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'Equations of the structural relaxation' in the Fi2-model are equations of motion de- 
termined solely by the vertices f i and V2 which play the role of the static structural infor- 
mation entering the MCT vertices V in eqn. (^). Bifurcations in the model lie on two lines 
parametrized by the exponent parameter A (Gotze 1984): type-B lines, where = 1 — A > 0, 
at vl = (2A — 1)/A^ and f 2 = 1/A^ with 1/2 < A < 1, and type-A transitions, where = 0, 
at = 1 and f 2 = A for < A < 1. Whereas in eqs. (|l|) to @ the microscopic transient is 
required in order for solutions to exist (Gotze 1991), after partial integration the short-time 
variation can be neglected, d'^^{t) + Jq dt' M'^^-{t - 1') df^it') < Vl?^{t), leading to: 

$(t) = mit) - ^ r dt' mit - t') $(t') , (6) 
dt Jo 

while the expression for the memory function is not changed: 

m{t) = vi<!>{t) + V2<^\t) . (7) 

Demanding the solutions to be of regular variation and that the Laplace transforms $(s) 
and rh{s) exist, enabling one to use Tauber and Abel theorems, see e. g. (Feller 1971), the 
equations of structural relaxation have to be completed by specifying the initial variation: 

$(t) ^ (t/t,)-'/^ for t ^ . (8) 

In this case, the convolution integral, which will be abbreviated as fQdt'm{t — t')$(t') = 
(m * $)(t), leads to (m * $)(t 0) 5(1/3, 2/3), where B{x,y) is Euler's Beta 

function. The Laplace transform of eqn. (P) for all s then is: 

{1 + s m{s) )^{s)= m{s) + V2t,B{l/3, 2/3) , (9) 

It is a straightforward exercise to verify that the last term in eqn. (H) is required for the 
existence of solutions with the requested properties (Voigtmann 1998), thus leading to eqn. 

Note, that the bifurcations to non-ergodic solutions cause $ and m to increase for small 
frequencies, suggesting to neglect (s + M^'^^-)Q^'^ <^ m{s) in eqn. (|^) in order to derive 

4 



equations for the structural relaxation alone. However, this procedure misses the non- 
analytic behaviour of {m * $)(t ^ 0) connected to the limit s —>■ oo. 

One can extend eqn. (|^) by a short-time series expansion, see appendix A, 

= it/Q-^/'' ( 1 + ^ c„ (t/Q^/^ ) for t < r . (10) 

n=l 

Figure 1 shows the numerically obtained radii of convergence, r, of eqn. (|10]) along the line 
of bifurcations in the Fi2-model. 

Figure 2 shows numeric solutions of eqs. (|^) to which exhibit slow relaxation processes 
upon approach of a bifurcation singularity at critical values of the coupling vertices. This 
explains the sensitive dependence of the dynamics on small regular changes in e, which 
measures the distance to the singularity. The algorithm described in (Fuchs et al. 1991) 
is well suited for the numerical integration. A type-B fold bifurcation at e = separates 
ergodic, $(t oo) = / = for e < 0, from non-ergodic dynamics, / > > for 
e > 0. Asymptotic expansions identify two divergent time scales connected with two power 
law relaxation windows (Gotze 1984). As demonstrated e. g. in (Franosch et al. 1998) 
the long-time solutions of the complete MOT equations of motion, eqs. (|l]) to (H), become 
independent of the microscopic transient. For long times they agree with the solutions of the 
equations of structural relaxation, eqs. (||) and (0). This is exemplified in figures 2 and 3 using 
two approximations for the regular microscopies in the Fi2-model: undamped oscillation, 
]\^reg. ^ pure relaxation, (s + M^^^' {s))Q~'^ — > Fq. It is apparent that oscillatory 

short-time dynamics masks the structural relaxation more strongly than relaxational. 

Assuming regular variation of the solutions of eqs. (|I]) to (^, at a bifurcation point, 
power law relaxation could be argued: 

$(t) = + /i (t/to)-" , for t ^ oo and wi,2 = , (11) 

where the exponent a follows from the exponent parameter A, which is a function of Vi2, 
via A = -8(1 — a, 1 — a) (1 — 2a) (Gotze 1984). From the equations of structural relaxation, 
eqs. (|^) and (^, and the short-time expansion eqn. ([T0|), one concludes that this asymptote 
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requires t > t^^r in general. The transient time to is obtained from matching eqn. (0) to 
eqn. (|1^); while the critical amplitude /i = (1 — /^) for the -Fi2-model. 

In the Fi2-model there exist two points where the short-time expansion can be found 
explicitly and solves the equations of structural relaxation, eqs. (^) to (||), for all times, i.e. 
r — > cx) in eqn. (]TD|); see appendix A and figure 1: 

$(t) = r + (t/t*)-^/' , for A = A, = I) = 0.684 ... and Vi^2 = <2 • (12) 

Both points are bifurcation points, as finite distances from the singularities lead to the 
existence of long-time relaxation rates as shown in appendix B, see also (Gotze 1984, Gotze 
1991). Furthermore, the initial condition, eqn. (j^), requires for the critical exponent a = 
a* = 1/3 which is satisfied for A = A,,,. Thus, eqn. (|12|) verifies power-law relaxation at a 
bifurcation singularity of equations describing the structural relaxation of a schematic MCT 
model. 

The existence of a solution of the integral equations for all times at a singular point also 
allows for these special cases to deduce a long-time solution for small distances from the 
bifurcation in an expansion which can be expected to have a finite radius of convergence, 

{t/to)\e\^/'<rp. 

m ^ {t/to)-'/' [ 1 + 5: A„ ( £ {t/tof' T], for t oo . (13) 

n=l 

The calculation sketched in appendix B, thus introduces the /5-time scale, = to 
a = 1/3 for A = A*, and agrees with the expansion of the so-called /^-correlator for small 
rescaled times t = t/t^; see (Gotze 1984). There it is also shown for general A, how to extend 
this solution to even longer times. Note that the expansion carried through in appendix B 
does not require a matching to the critical decay. 

§4. Extensions to MCT equations for simple liquids 

The extension of the analysis of the bifurcation dynamics of schematic models to the full 
MCT equations for simple liquids, eqs. (0) to (^), has been achieved in (Gotze 1985) and 
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rests upon the central-manifold concept. The generalization of the equations of structural 
relaxation to coupled, wave vector dependent equations is trivial; eqn. (||) attains g-indices 
and eqn. (|^) is replaced by eqn. (^). For quadratic polynomials in eqn. (^), the power law 
of the initial condition is not changed: 

<^g{t) ^ Xg r^'^ for t ; X, > . (14) 

It is surprising that a vector Xq has to be specified in order to render the solutions unique; 
recall that owing to the factorization property of the critical decay, ^q{t) = fg + hgit/to)^"- 
(Gotze 1985), the long-time dynamics of the structural relaxation becomes unique by spec- 
ifying only the one matching time to- 

From simple generalizations of the calculations presented in appendix A it is obvious that 
the short-time expansion, eqn. (|l3), can be obtained in the wave vector dependent case also. 

§5. Conclusions 

Definitions of equations for structural relaxation not restricted to long times are given 
within the MCT. These equations, eqs. (|) and (||), together with the initial conditions, eqn. 
(0), determine smooth but sensitively varying intermediate scattering functions $g(t) from 
the equilibrium structure of the liquid; the static structure factor Sq is the only quantity 
appearing in the equations of motion. The possibility to specify an arbitrary initial non- 
negative amplitude, Xq, is worth stressing and also holds for the discrete dynamics model 
of (Franosch et al. 1998). The short time divergence of the solutions of the structural 
relaxation equations renders them unphysical for too short times; they violate the imposed 
normalization of the $g(t). This stresses that the physical mechanism (the so-called 'cage 
effect') causing structural relaxation is relevant for long-time dynamics only. Any complete 
description of the density fluctuations requires a transient to regular microscopic short-time 
dynamics. Nevertheless, extending the structural relaxation equations to all times on the 
one hand allows the deflnition of a lower limit of validity of this long time description — 
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the time where the solutions increase above unity could be chosen — and on the other 
hand presents a conceptual simplification as some transient effects are eliminated. Thus for 
example, analytical proofs of power-law relaxation at bifurcation points become possible for 
special cases; recall that at regular points of the full MCT equations a longest relaxation time 
exists. Also, the smaller of the two divergent time scales of MCT could be identified from a 
perturbation expansion. Moreover, the aspect that MCT equations of motion with different 
models for the short-time microscopic motion lead to the identical long-time dynamics is 
again observed. 
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Appendix A 

Entering the ansatz = + (pit) into eqs. (|[) and (0) leads to 

a = a (Pit) - A (i)\t) + (1/A' + a) (0 * 0)(t) + (02 * 0)(t) (Al) 

where a dot denotes a time-derivative and the following abbreviations are used: A = A/ A', 
A' = A/(l - r), a = (1 - r)V/A, a = (1 - ni6v,r + 5v2r'), A = (A' + (1 - n'6v2), 
1/A' + a = [1 + (1 - f^)\6vi + 26v2f^)]/A, and a = (1 - /^)(1 + a - (1 + a) /A) which are 
obtained from Vi^2 = "^12 + ^'^1,2, where the limit of V2 approaching zero has been excluded. 
In the vicinity of a type-B transition. A' = 1, A = A. Close to a type-A transition A' = A, 
and thus A = 1. From the short-time expansion, eqn. ([To|), one obtains a recursion relation 
for Cn, where Cq = 1, Ci = Ci — and c„ = c„ else, and 6n,m is a Kronecker symbol. 

Cn+l = fr,T^(n+2^2\ , j(n+3^T\ { ^[^^ ^^)(1/A' + (x)] Cn-mCm + CT6n,2 " 

3 ' 3^ ~'~ V 3 ) 3JJ m=0 



n m 



;i - M[^C„.i + ^ E (E C^-kCk + l)Cm)~Cn+r^m] } (A2) 

m=l fc=0 



At two singular points, where & = a = a = and A = A* = |)/3 = 0.684 . . ., all c„ 
vanish. This proves eqn. (|l^). Approaching the Fi-model at (fi,f2) = ("^1,0) (Gotze 1991), 
the radius of convergence of eqn. (p!0| ) shrinks as for a m-th order polynomial mode coupling 
functional a similar expansion with (t/t*)^^/^™^^^ for t ^ applies (Voigtmann 

1998). 

Appendix B 

Close to the type-B transition at A = A*, the eqn. (|A1|) can be expanded around the 
solution 0o(t) = t~^^^ at 5vi^2 = 0. Using 0(t) = 0o(t) + 0i(t), specifying the path such that 
5" = for simplicity and integrating this leads to: 

ds Kit, s) 0i(s) =at- lat^/^ - ((a + 0i + 20o0i) * 0i)(t) + ((A* - <po - 0i) * 0?)(t) , 

(Bl) 

where the Volterra kernel is of generalized Abel form: 

K{t, s) = 2{t- s)-'/' - 2K s^'/' + (t - s)-2/3 + 2 (t - s)-i/3 . (B2) 

From the short-time expansion eqn. (|10|), (f)i{t ^ 0) oc t^^^ can be deduced in this case, ruling 
out that Kernel K has eigenvalue zero. The expansion 0i(t) = Y,n=i d^^K^) thus leads to 
a recursive set of solvable linear integral equations for the g^"'\ If their long-time behaviour 
is searched for, only the first two terms in K need to be considered and by induction one 
shows that c/(")(t) ^ An t(2«-i)/3 for t oo, where the An obey: Ai = 1/(2(5(|, |) - A,)) 
and 

1 n— 1 

= [ E ( ^) - A. ^ ) 2n±2) _ A, 2^ ) . (B3) 

Replacing a with e, this coincides with eqn. (^) and the short-time expansion of the /3- 
correlator for A = A* (Gotze 1984). 
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Figure captions 

Figure 1: Radius of convergence r of eqn. (p!0|) numerically obtained using the Cauchy-Hadamard 
criterion and eqn. (|A2| ) along the line of bifurcations in the Fi2-model. The line 
parametrized by A is given by V2 = and vl = (2A — 1)/A^ for type-B (/^ > 0), 
and by f 2 = A and Vi = 1 for type-A (/^ = 0) transitions. The vertical dashed lines 
indicate A* where eqn. (|TU|) with Cn>i = holds for all times; i.e. r 00. The dot- 
dashed line describes the linear decrease, r oc A, see appendix A, upon approaching 
the Fi-model at (^^,^2) — (l^O)- 

Figure 2: Correlators ^(t) obtained from eqs. (^) and for exponent parameter A = 0.70. 

The bold line marked by a c corresponds to the bifurcation point and the others are 
calculated for t>i_2 = 2(1 + = ±1/4", = 0, 1, . . . as labelled. The time scale 
to from eqn. (|l^) is indicated and used in the inset to compare the curve at v'^ to one 
crossing over to oscillations at short times (dashed line calculated for the Fi2-model 
from eqs. (|l]) to (Q) with M^^^' = 0) and to another with overdamped short-time 
dynamics {{s + M'^'^^- {s))Q^^ —>■ Fq, long dashes). 

Figure 3: Comparison of correlators at bifurcation points of the Fi2-model for different micro- 
scopic transients and shifted using are plotted versus rescaled time, with to from 
eqn. (p!T|). Solutions for exponent parameters A* = 0.684 and A± = A* ± 0.15 as la- 
belled are shown obtained from the equations of structural relaxation, eqs. and 
(solid lines), and from eqs. (|l]) to (^) with undamped oscillatory (iVf'^s- = g, short 
dashes; lines terminated for t^lO^to) and with overdamped ((s + M'^'^^- {s))Q~^ Fq, 
long dashes) short-time dynamics. Dashed-dotted curves indicate the critical decay 
law, (to/t)"' with a+ = 0.26 and a_ = 0.39 corresponding to the A±. For A*, where 
a^, = 1/3, the structural relaxation follows the critical decay for all times. Black circles 
(diamonds) indicate where the structural dynamics (critical decay) curves deviate by 
20 % in horizontal direction from the solutions with overdamped microscopies. 
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